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ABSTRACT 


The Dynamic Model of Modern Military Conflict developed 
by Dr. Paul Moose, Naval Postgraduate School, is described 
by its system of differential equations followed by an inves- 
tigation of its behavior. This investigation is predicated 
by an analysis of the model's stability about equilibrium 
using a method attributed to the study of ecosystems. The 
basis for this analysis is the formulation and subsequent 
evaluation of a community matrix termed the conflict matrix 
men this investigation. 

Specific investigation of the model's behavior is observed 
varying the Command, Control, Communications and Intelligence 
ie o) enhancements of one combat force against that of an 
opposing force while one force utilizes misinformation and 
deception as a counter 3 tactic and the other uses physical 
destruction as a counter C? tactic. 

Potential military use of this model as an analytical 
tool for playing "what if?" type wargames is envisioned 
following further research and study of the effects on model 
beNavior by varying other parameters not specifically addressed 


moeeciis initial investigation. 
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it Ln obUuc Tl LON 


A military conflict has a variety of elements that 
potentially influence its eventual outcome. Commanders of 
the various echelons within the combat environment determined 
mea particular military conflict are faced with decisions 
that require the fusion of information from countless sources. 
These sources may be highly sophisticated sensors, relaying 
timely digital information concerning enemy force concentra- 
Pons Or moyements, aircraft or missile tracks both friendly 
and nonfriendly, or perhaps intercepted communications signals. 
These sources could also be the professional schools, exercise 
participation or intelligence documents that have molded the 
commander's decision process over his military career whether 
it be extensive or minimal. Regardless of the case, decisions 
are made with some reliance on or rejection of prior knowledge. 
The assimulation of this prior knowledge, whether relative to 
the situation or not can be a monumental task in itself. 

The sophistication of computer driven information systems 
and a presumed need for rapid decision making either by the 
man aS a commander or the machine as his extension has caused 
a proliferation of model development in recent years [Ref. l]. 
These models attempt to reflect the variety of elements that 
potentially influence the outcomes of military engagements, 


manipulate them according to some predetermined scheme or 
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algorithm and finally output some prescription or calculated 
result. The output of such models, together with the parti- 
cular input elements and scheme of manipulation offer 
limitless applications either for direct command decision 
making, training or running countless “what if" exercises 
relative to known or suspected enemy and friendly combat 


mimreornmation or characteristics. 


pee: URPOSE 

The primary purpose of this thesis is to initiate inves- 
tigation of one such model that exhibits the potential for 
Mere In a variety of areas, but in particular in looking at 
the “what if" questions associated with command, control and 
communications (C3) and counter c3, The specific model under 
mivestigation is the Dynamic Model of Modern Military Conflict 
proposed by Dr. Paul Moose, Academic Associate for the Joint 


Monmand, Control and Communications curricula at the Naval 


Postgraduate School [Ref. 2]. 


eee SCOPE 

. miemoasic Gescription of the model is given in Chapter 2. 
This description is in the form of a system of differential 
equations similar to that which describes Lanchester-type 
attrition models, however it includes terms relative to c3 
and counter c? tO: esposiung forces, in this respect, informa- 


tion possessed by opposing forces is elevated to an equivalent 


level with the forces themselves. Furthermore, the description 


1 





includes a discussion on how the equations of this model are 
also similar to evolution equations which descrihe multi- 
species ecosystems. Such systems include natural death (loss) 
rates and replenishment of species within a population thus 
leading to a mathematical discussion of this model's behavior 
relative to an established equilibrium. 

Following this description, Chapter 3 looks at the 
Beability of the model in terms of a “community matrix", 
another reference to ecosystems which incorporates those 
parameters which tend to increase or reduce species within 
peeooDuUlation Or community. This “community matrix" is 
renamed the "conflict matrix" for this model and explored 
Dy asSigning a certain set of asymmetrical values to the 
elements of the matrix relative to hypothetical opposing 
forces in a military conflict. Specifically the investiga- 
tion undertaken in this chapter, demonstrates conditions 
where opposing forces stabilize at a set equilibrium point 
or diverge from this point in unstable fashion. 

Chapter 4 looks at more specific behaviors relative to 
Bette stability and instability demonstrated in Chapter 3. 

The identical initial conditions used in Chapter 3 are 
again used in this chapter and represent arbitrary units 
rather than specific calculated values. These values were 
chosen purely to demonstrate the general behavior of the 


moecel rather than replicate any specific military conflict. 


IL a 





Pindilivy ckrects On equilibrium relative to the initial 
conditions used during stability and behavior investigations 
eee Explored in Chapter 5. Specifically, the initial condi- 
tions are aimed at isolating the effects of varying the 
command, control, communications and intelligence (C?r) 
enhancements of one force against an opposing force's cot 
enhancements. Utilizing these results, further demonstration 
of the model behavior is then achieved by varying resource 
and information variables associated with opposing military 
merces. 

In essence, this thesis is concerned with an initial 
demonstration that the Dynamic Model of Modern Military 
Conflict can perform as might be expected rather than prove 
that it works in all cases. This approach opens several 
avenues to continued research leading to a more complete 
demonstration of its usefulness in evaluating the influence 
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=e C L and counter ee iieMedern Military conElices 
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It. A DYNAMIC MODEL OF MODERN MILITARY CONFLICT 


fee DEFINITION 
The Dynamic Model of Modern Military Conflict incorporates 


$y Wikehiisds1C eater tion COecEfticients such as 


meme effects of C 
those found in Lanchester-type models. Specifically the 
effects of the above are described in terms of information 
(IL and ry) and general purpose resources OS and M) fOr 
opposing xX and Y combat forces [Ref. 3]. The system of 
differential equations which characterize this model are 
gaven below followed by sections in this chapter which 


briefly touch on the background of the model and provide 


specific definitions of the terms involved. 


GPiemy 2.) 
d(I,)/dt = I,[-A, - _— Ye) eee Le 
(Hames 2 <2) 
d(M.) Vat := M.[-B,, - Sxyty - By - aM At - Fel y = De aa RAM 
(Eqn 2.3) 
aI Gua | |=} =f [= M +C JT M + 6 
(T)/ yh Ay 7 yal 7 Vath + Stell QTE 
(Eqn 2.4) 
oem )/dt =M [3 -6 I = M]=-$dM I-dMI=-bM M +RM 
y) /at yl Vann ace rae yy yey yx ye x “yx yex -‘y ye 
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B. BACKGROUND 


This system of four first order non-linear differential 
a(S) 
ate 


equations in this form are known as evolution equations. 


equations is of the form A set of 


mia A 


232 


S and Q are 4xl column vectors and there are four func- 
tional relationships which are considered to be at most 


quadratic. Expanding the loss function term gives the 


moLLOwing: 
Geen 2. 5) 
4 4 
- F, (Ss) = 7830 295585] oy a,.8;, i= 1,2,3,4 
ie. 
a F. = - quadrative term - linear term (Sqn 2.6) 


This system as it turns out is a very generalized form of the 
Lotka-Volterra equations for multi-species ecosystems [Ref. 


4: 37-38]. More specifically it takes on the form: 


CE Gine2. 7) 





4 4 
ale ae 2 es 
rae) aie ate p eS. harris efOa 2 = 91,2534 


fe EAPLANATION OF TERMS 


In general the greek letter coefficients are quadratic, 
and the others are linear. Following is the definition of 


mee terms that characterize the model: 


I. ~PcOEMAtloOneor LOorce X 
= - (Same as above at equilibrium) 
i General purpose resources for force X 


iS 





M - (Same as ahove at equilibrium) 


xe 

A, - Natural death (loss) rate for OT resources 
Gt EOrce 4 

Be - Natural death (loss) rate for general purpose 
resources of force xX 

B rx - Lanchester "area fire" coefficient, X on Y 

ve Lanchester "aimed fire” coefficient, X on Y 

 - Fixed replenishment rate for ea resources of 
hORee... 

i - Fixed replenishment rate for general purpose 
resources of force X 

a - Replenishment from general purpose resources 
Pomme o© PprEeaduction, £once X 

= - Diminishment of X force for crt production 

= Coun cer os by force X on force Y by misinformation 

Y and deception from C3I resources 

Regs: = Counter ce DPymetonGe.. On force Y by physical 

Y destruction from general purpose resources 

6 - Gen (quadratic) force effectiveness enhancements, 

yx 
X on Y 

a - cy (linear) force effectiveness enhancements, 


ve SPO} ae 4 
The remaining coefficients with y subscripts or xy subscripts 


Meerter to the Y force or Y on X activities. 


Pee GAPECTED MODEL BEHAVIOR 

The actual behavior of the model is dependent initially 
mmeene equilibrium for the system. One equilibrium point is 
determined by solving the following linear set of equations 


ee and M 


ror Tyee! M ye ye’ 


AS 
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(Eqn 2.8) 


: “304 “xy ~Txy rxe] | x *x 
=o Q = Std ay) = (8B, tb.) Re ee| eR 
=O - 0 Cc i A = 0 
yx Vyx yy ye ley | a 
-(S6 +d_) -(B _+b_) -d 0 M B R 
yx YX yX “yx yy ye y y 


Once equilibrium is established the effects of varying 


ase Mos qe a 


can be illustrated by solving the system of differential 


certain coefficients on dependent variables I 


equations over the independent variable time. 





Mees bee ety ANALYSIS 


eee Der INITION 

Stability, as it relates to a deterministic model or 
deterministic system of equations, is more specifically 
defined as "neighborhood stability" or "stability in the 
mmeimity OL equilibirum”. For population models with 
environmental parameters all well defined constants, 
interest centers on the community equilibrium where all 
species’ populations have time-independent values, i.e. 
all net growth (death) rates are zero. Such an equilibrium 
may be calied stable if, when the populations are perturbed, 
they in time return to their equilibrium values; the return 
may be achieved either as damped oscillations or monotonically. 
Conversely, if such a disturbance tends to amplify itself, the 
system is called unstable; such instability may appear as an 
oscillatory or as a monotonic growth in the disturbance. 

This definition is sufficient for linearized systems, however, 
ag may be misleading for non-linear systems when full global 
meaobility of a system is of concern [Ref. 4:15]. 

Full global stability is said to be characterized by 
neighborhood analysis if the existence of a Lyapunov function 
can be determined [Ref. 4:15]. The existence of such a 
function is not explored in this thesis, however, regardless 


of whether such a function exists or not, neighborhood analysis 


18 





is within the scope of this thesis and founded in the 


generation of a community matrix. 


Be COMMUNITY MATRIX 

The community matrix is said to both summarize the system 
(its elements being determined by the interactions between and 
within species near equilibrium) and sets the neighborhood 
stability by the sign of its eigenvalues (all negative real 
parts for all eigenvalues imply stability). Specifically it 
is generated by the following mathematical derivation [Ref. 4: 
I-22}. 


Since the model is of the form 





= -F(S) + Q 


its multispecies population dynamics are given by the set 


Be LOur equations: 


(Eqn 3.1) 
a(S, (t)) | 
—ge = FYZCS) (4), So(t), S3(t), S,(t)), i = 1,2,3,4 


“The growth rate of the acs Species at time t 1s given by 
the nonlinear function F,; of all relevant interacting 
* 
Beepulations. Equilibrium populations, S; PD Peeleuew £rom 
the 4 algebraic equations obtained by setting all growth 
Maees at zero: 
(sen. 322) 
* * 


Ike, 





Her each population, let 


* 
S, (t) = S. “- s. (t) (Hane 32.3) 
where S; measures the initially small perturbations to 
the ee population. Expanding each of the basic equations 


around this equilibrium in a Taylor series and discarding all 


second order and higher terms in the population perturbations 


Ss, a linearized approximation is obtained: 
d s. (t) 4 
a 5 ais 2. te) (Eqn 3.4) 
gl 


Equivalently, in matrix notation 


GaeGE) 
Cue 


where X is the 4xl column Matrix of the Ss and A 1s 


= A X(t) VSM es 45) 


the four by four "community matrix" whose elements ace 
describe the effect of species Jj upon species 1 near 
equilibrium. The elements S43 depend both on the details 
of the original equations and on the values of the equilibrium 
mepulation through 

Bhd e Sg 


= a (ae) (Eqn 3.6) 


3) 


The partial derivatives, denoting the derivatives of F 
keeping all populations except Ss. constant, are then 
evaluated with all populations at their equilibrium values. 
The corresponding community matrix for the Dynamic Model of 
@emerhn Military Conflict is called the “conflict matrix" 


and its elements are given in the following matrix: 
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eee a  Y ) @ —o -Y 


xy x xy xy 

= Sheen! TS ge ce a dace 
“A Yyx APY) Sy 

“(6d ) -(8,¢B ~i,,, (B45 +8. 


ee METHOD OF ANALYSIS 

To investigate the stability of the model, eigenvalues 
are computed relative to the conflict matrix given above and 
a set of initial values for the coefficients that comprise 
the elements of the matrix. Two cases are investigated: 
one in which element Coxe is held constant and element dey 
1S varied across a range of values and another where Gay 
is held constant and Cs 1S varied across a range of values. 


The particular method of computing eigenvalues is the IMSL 


routine EIGRF [Ref. 5]. 


eee LNITLIAL CONDITIONS 

The initial conditions for investigating the two cases 
are given in Table I. These initial conditions were chosen 
to illustrate some of the potential behavior that this model 
may exhibit. The variable parameters are err enhancements 
for both xX and Y forces. Furthermore the model is 
asymmetrical; xX is using physical destruction by forces as 


S 


a counter C”~ technique and Y is using misinformation and 


deception as a counter om technique. 


Juak 





ALANS OE) gaa E 


INITIAL CONDITIONS FOR STABILITY INVESTIGATION 


A =A =B8B . =8 . = Q.5 
x yY x yY 
oe = SK = oon = oy = Q.0 
Py = Pox = 1.0 
6 = 6 = 1.0 
XY yx 
b =) 5 = 0.0 
XY Vx 
eet = 1.00 (misinformation and deception) 
Ox = 0.0 
Y xy = 0.0 
vx = 1.0 (physical destruction) 
Gay = variable cx enhancements (Y on X) 
d = variable ctx enhancements (X on Y) 


vex 
EGuieanium @o taken GO be the vector 1 , i.e. 


1 = M = I = M = 1 
xe xe ye ye 
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SisctteultioOnm OL che initial parameter conditions into 


the conflict matrix gives the following matrix: 


Ba Q =-1.0 0 
0 =i -d = 0 
yx 
0 -1.0 -1.5 0 
“dose il 210) 0 =Le5 


The particular cases investigated then are: 
1_ Aas | held constant at l andd varied from 0 to 2 
yx ory, 
by 7Om> . 
Zee OC. held constant at 1 andd VaGrecd trem 0 EO 2 
AY We 
Dy 20 20'5.. 


Peek oULTS AND DISCUSSION 

The eigenvalues generated by case l are presented in 
Table II. It can be seen that for values of fay greater 
than 0.80 a positive real eigenvalue appears. This is 
mi@mcative Of instability. Similarly all real parts of 
eigenvalues for dey less than or equal 0.80 are negative 
megs indicating stability. 

The eigenvalues cenerated by case 2 are presented in 
Meese III. For values of d G@mcarereenanieOl. oo) InStabilicy 


YA 


mempindicated while stability is indicated for values of ax 


Mees than or equal to 0.55. 
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TABLE If 


EIGENVALUES FOR VARIABLE en ENHANCEMENTS (Y ON 4&) 


EIGEN 3 EIGEN 4 


EIGEN 2 


EIGEN 1 


DXY 
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The results of the ahove two cases are summarized 
Prapnically in Pigures 1, 2, and 3. Figure 1 depicts the 
generation of eigenvalues with. positive real parts 
eenstability) beyond 0.5 and 0.8 for dx and ae 
respectively, and negative real parts elsewhere (stability). 
Figure 2 represents the other real eigenvalue which is 
always negative thus not contributing to instability. Figure 
3 represents the positive imaginary part of the remaining two 
eigenvalues which turn out to be complex conjugate pairs with 
real part equal to -1.5. Existence of these imaginary parts 
indicate that the system does in fact oscillate while dampen- 
mag Out. 

The eigenvalues generated by these two cases lead to 
more specific investigation of stability and instability. 

The eigenvalues have located the regions of stability and 
instability so now by using the appropriate values of day 
and tho the behavior of the model can be observed relative 


to perturbations of the dependent variables over time for 


stable and unstable cases. 
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Figure 2. Eigenvalues with all negative real 


parts (imaginary part = 0). 
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Figure 3. Positive imaginary parts of complex 


conjugate pairs (real part = -1.5). 
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ee Sewer C SrABTLITY ANALYSIS 


A. METHOD OF ANALYSIS 

More specific illustrations of stability in the Dynamic 
Model of Modern Military Conflict are achieved by perturbing 
the dependent variables and observing the behavior of the 
model relative to time. The initial values used in the 
stability analysis of Chapter 3 are substituted into the 
differential equations which characterize the model. These 
differential equations turn out to be the following for the 


parameter values of Table I. 


ole 
x = = 
Soa eS t, (I. a (Gs) (Eqn 4.1) 


dM 


x = — = 
= (| laa oer M. (M,, + 0.5) + ae (Eqn 4.2) 


ee 


a = 
it eS I. Oh ee s5) (Scie 2) 


dM 
=1.5 -d I. - M ie OS etal (Eqn 4.4) 


she WPSae.s yx 
-The Interactive Ordinary Differential Equations package 

(IODE) was employed using these differential equations and 
representative values for do and oh (stable and 
unstable values) [Ref. 6]. Specifically, plots were generated 


to graphically represent two behaviors: 


1. Dependent variables against the independent variable 
meni, Lea, M against time. 
x x y Y 


eee nase Plots Of Opposing forces i.e., ne against M,, 
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Be ola opr LiTyY WITH dy AT Q.4 

Table IV shows the initial values for the dependent 
variables, time duration and increment for plotting for 
each of fourteen trial runs of this model. Both ay and 
dose were held constant at Q.4 and 1.9 respectively. In 
addition, the final column in Table IV represents the time 
at which all dependent variables damped out to or very nearly 
to the estahlished equilibrium of l. 

Graphically, Figure 4 depicts the general plot where 
only one dependent variable was perturbed (trials 1 through 
wee eilogures 5 through 8 further illustrate phase plots for 
these same trials of My vs Me PMustwatcing return to 
equilibrium over time. 

mvomeener trials ace illustrated in Figures 9 through l2. 
Specifically Figures 9 and 10 represent perturbations of both 
I, and (trial 6). The interesting phase plot in Figure 10 
@s indicative of the oscillatory nature of both M and Me 
relative to these perturbations. Figure ll depicts a triple 
Perturbation involving Lie I, and Bho with the resulting 
SOscillatory nature of M, and M., aroun ebugure 12 (Trial 13). 

Although the above illustrations are not all inclusive, it 
Should be noted that all dependent variables did in fact damp 


out over the set time period and phase plots of M,. and MS 


revealed return to equilibrium in all trials. 
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TABLE IV 


INITIAL CONDITIONS WITH 


he @.4 AND qx 
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See oLABLILITY WITH qx Py (Oe EO 

Similar to Table IV, Table V lists the initial parameters 
for investigation of the model with qx aemO.10- and dey 
meee.00. Graphically, for a single perturbation, this form 
of the model produces similar results as Rezore The excep- 
mien 1S given in Figure 13 where. only I. 1s perturbed 
Se~ial 3). M, and M can be seen to return to equilibrium 
without oscillation as had been indicated in the previous 
investigation. Again all trials exhibit the dependent 
variables damping out in the time allowed and since no 


Significantly differently graphical forms were produced they 


mee not imcluded in this section. 


Pei NoTABILITY WITH Qoey AT 1.40 

Table VI lists the initial parameters for the unstable 
case where doy = 1.40 and Gos is held at 1.00. The 
additional column on this table represents which force term 
M. or My 1S increasing without bound while the other term 
1s decreasing. It should be noted that some variables were 
given initial values of 1.001 to get comparative plots over 
time, 1.e. using 1.25 as an initial parameter in these cases 
caused either M. or MS to increase or decrease severely 
immediately. 

Graphically three general forms of instability were 
demonstrated where only one dependent variable was perturbed. 


fiiese are presented in Figures 14, 15, and 16 (trials l 
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TABLE VI 


INITIAL CONDITIONS WITH 


d AT 1.4 AND d AT 1.0 (UNSTABLE) 
XY yx 
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EProugi@ea). Ine corresponding phase plots are given in 
Peaures 17, is, 19, and 20. Specific values for double 
and triple perturbation did not demonstrate any unusual 
graphic forms of instability and generally followed those 


graphic representations displayed for single perturbations. 


Pee LNSTABILITY WITH Gs AT 1.40 

Mable VII gives the initial values for trial runs holding 
doy eee OO and oe at 1.40. Similar graphic forms were 
produced in this investigation as those produced in the 
preceding investigation of instability. Two exceptions are 
given in Figures 21 and 22 (trial 7) and Figures 23 and 24 
(trial 10). It also should be noted that in this case. 
perturbing M,. and IY together causes M, to increase. 


In the previous case perturbing these variables together 


drove a upward, an interesting result. 


F. GENERAL OBSERVATIONS 

The cases investigated in sections B and C of this 
chapter clearly indicate that the Dynamic Model of Modern 
eiitary Conflict does exhibit stability as predicted in 
Chapter 3. Furthermore, it also may behave in an unstable 
manner as indicated in sections D and E above. A closer 
look at the unstable cases also reveals bifurcation relative 
to the values for day and those! tical values for the 


Gepencent variables and the resulting increasing term. 
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A simple display of this bifurcation is seen when only 
one element (dependent variable) Is perturbed as in trials 
1 through 4 in Table VI and VII. Relative to Table VI, doy 
is set at 1.40 and qx Spee OG. fais Gives ¥ an 
advantage over X in C3 effectiveness, however, it can be 
seen that in trials l and 2 of this table, X can overwhelm 
fe (ML increasing, My, decreasing) if given an advantage in 
either information (IL) Sr rOrces (M.). Similarly, relative 
to Table VII with dy set at 1.00 and See ace 1.402 ¥ 
can be seen to overwhelm xX in trials 3 and 4. Evidences of 
bifurcation are a bit more complicated in the remaining trials 
of these tables. 

mimtablo Yio trials 5, 9, ll, and 12 exhibit this 
bifurcation where y has an advantage in a but is over- 
Pmemmea by XxX. Similarly in Table VII, trials 6, 7, 8, 13 
and 14 show xX being overwhelmed by Y, even though it has 


a c? 


advantage. Trail 10 is not a bifurcation example, but 
as mentioned in section E, does produce an interesting result. 
Beecifically, from Table VI, a slight advantage in MS 
-and L., results in Y overwhelming X when Y has a c 
advantage. Whereas from Table VII, the slight advantage in 
x 


M. and I, results in xX overwhelming Y when X has the 


ce advantage. This seems consistent with common sense. 


ao 





V. EQUILIBRIUM ANALYSIS 


All the preceding analysis of stability was predicated 
upon a constant equilibrium vector of 1 for ease of calcula- 
tion. Another question that should be answered at this point 
is "what happens to equilibrium values if, for instance, when 
replenish rates (ORY QR) are kept constant and a and/or 
ope is allowed tO vary?" To answer this question the equi- 
librium equations and initial analysis conditions are brought 


back into the analysis. 


A. METHOD OF ANALYSIS 

First, the original equilibrium equations are set to 0 and 
solved for the replenishment rates. This result is then used 
to solve the same set of equations for equilibrium as a func- 


mien of a and a 
cy yx 


PeeeesOUTLLIBRIUM RELATION TO 4d AND d 
XY yx 


Each equilibrium equation is set to 0 and solved for 


replenishment rates. 


Cx te 7 tyytye 7 Yey ve Gree oe Tk 
Oeil) = (1) (1) - O(1) - 0.5 + Q., =) 
oF = 1.5 (ig aero.) 
Be xe 7 a . (Bey : eee 7 Be Re 
O(1) - (0 + Gy) IL —l + 0) l1-=+0.5 + R. = 
Ry = os + 1.5 (Gan. 5n..<) 
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Bo ee + Cyy Lh 0 
(aaa —~ ttl) + ) Oil) - 0.540, = 0 
Q, = 1.5 CEG = 5. 3) 
“(6 + Gn) te - (Bo a Do! Mee - doytye - BY - RY = 0 
Sede CO ub - ee, ) 1 - 0 (1) - 0.5 + R_ = 0 
( ye) ( ) (1) : 
Be = dos + 1.5 Geena a4 ) 


Utilizing these results the same equations are solved 
again for the equilibrium points starting with some arbitrary 
—. . . ; 
initial values for dey and Gone! day and antares respectively 
to fix the replenishment rates. 


CM -a I - y M -~ A +Q = 0 

xx xe xy ye xy ye <x x 

0 ee) = (2) oe! - Os) = 0.5 + om = 0 
lve = Q., - 0.5 


Using the results of equation 5.1 


re = ] ag ne SiS) 
i an ~ On + Gey) Ive ~ (Boy + Oy ve - By + Ra = 0 
mee) t_. = 10) se dey Ive = Beer + 0) Moe - 0.5 + Bae = 0 
Moe = ~Gevive <0 Sear Ry 


Using the result of equations 5.2 and 5.5 and substituting 


| 
dy + 1.5 for Ry. 


M = -d =O Peoee tcl - 94 1. 5) 
ye XY XY 


M 


T 
= 
O° 


a 1 
ye ‘ele oe ) (Bqme> <6) 
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“txt xe ~ Vox ee Coy tye > AY ‘ 2, = 8 
= (O) TL. - (1) My. ~ (0) M. - q.5 FO, = 0Q 

Mee = Q,, = 3 
Using the results of equation 5.3 

Mo = | CEG y ) 
-(6 +d )I - (8 +b )M -d I =- 8B + R. = 0 

yx yx xe ae yx xe yy ye y y 

-(Q + don) i xe - (1 + Q J - (0) I. - Do ae Is = Q 


Using the results of equation 5.7 
or a 
y 
xe cl 
yx 


Using the results of equation 5.4 and substituting 


aoe t 1.5 for R 
y 


yx 
Cue ee de 
i = ee ee 
xe d 
yx 
Adding d -d to the numerator 
WK yx 
(d - d'_) 
iE = 1 — os) eee 
xe d 
yx 
ASSIMILATION 


It suffices to say then tnat Loe and Moe are not 


effected by changes in d or d=. On the other hand 


XY yx 
varies in the opposite direction to deviations from 


mesmnominal value. te varies in the opposite direction to 


deviations from its nominal value normalized by ds. 


ye 
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More specifically, equation 5.6 indicates that force Y 
would require fewer general purpose resources to hold force 
X to 1 provided, force Y increased its c? enhancement 
beyond its initial value. In a similar fashion equation 5.8 
indicates force X would not require as much information to 
hold force Y to 1 provided force xX increased its ce 
enhancement beyond its intial value. It should be noted 
that these equations may not be true in general, but rather 


may only be true for the specific initial values used in 


this analysis. 
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VI. SUMMARY AND CONCLUSIONS 


This thesis has described the initial investigation of 
the Dynamic Model of Modern Military Conflict. Following a 
brief description of the model, investigation of stability 
was undertaken using a method of analysis previously prescribed 
for the study of deterministic systems such as ecosystems. 
Through this investigation it has been demonstrated that the 
model does in fact behave as one might expect, given a set of 
somewhat reduced initial parameters. In this respect, this 
research endeavor did not attempt to fully prove or disprove 
the model but rather its purpose was to permit observation of 


the model behavior as an initial step toward further study. 


A. GENERAL OBSERVATIONS 

Tt can be seen that the Dynamic Model of Modern Military 
Conflict, given a set of asymmetrical initial parameters 
exhibits behavior that corresponds to a status quo (stability) 
or a dominance of one force, combat force if you will, over 
another (instability). The significance of this observation 
lies in the plausible use of this model in predicting outcomes 
@eettary Conflict. Particularly, the incorporation of 


3 and c3 


terms relating to counter C enhancements along with 
more traditional terms such as loss rates, replenishment 


rates and Lanchester coefficients, adds a contemporary flavor 
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to this model. The actual employment of this model in 
playing "what if?" type wargames however cannot be undertaken 
as yet without additional research including the remaining 


terms of this model. 


oe LMPLICATIONS FOR PUTURE RESEARCH 

Specifically, the following areas for additional research 
are suggested based on the initial investigation described in 
this thesis: 


1. Location of other equilibrium points besides the 
vector 1 utilized herein. 


2. Determination of an appropriate Lyapunov function 
Suitable for predicting global stability. 


See CCGQlCeMon OL Difurcation through an analytical 
approach. 


4. Further investigation of model bDehavior relative to 
variations of other initial parameters (i.e. counter C3 
Secrticients). 
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